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Abstract 

We consider generic models of quintessence and we investigate the influence of massive neutrino 
matter with field-dependent masses on the matter power spectrum. In case of minimally cou¬ 
pled neutrino matter, we examine the effect in tracker models with inverse power-law and double 
exponential potentials. We present detailed investigations for the scaling field with a steep expo¬ 
nential potential, non-minimally coupled to massive neutrino matter, and we derive constraints on 
field-dependent neutrino masses from the observational data. 
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I. INTRODUCTION 


It is remarkable that both early and late time completions of the standard model of the 
universe include phases of accelerated expansion. Inflation [1-3] not only successfully ad¬ 
dresses early time inconsistencies of the hot big bang, such as the flatness and horizon prob¬ 
lems, but it also provides the mechanism of primordial perturbations needed for structure 
formation in the universe. On the other hand, the age crisis [4] in the standard cosmological 
model again asks for a late-time cosmic acceleration for its resolution in the standard lore. 
This phenomenon was directly conhrmed by supernovae observations in 1998 [5, 6], and it 
was indirectly supported by other probes thereafter [7-9]. Such a conhrmation for inflation 
is still awaited. 

The late-time cosmic acceleration can be caused by the cosmological constant [10], with 
an energy scale of 0(10“^) eV, or by a slowly-rolling scalar held (quintessence) with the 
mass of the order of Hq ~ 10“^^ eV [11-14]. Both options are plagued with deep theoreti¬ 
cal problems. Concerning the cosmological constant, it is an “unnatural” parameter [15] of 
Einstein theory, and its small numerical value at the classical level gets destabilized by quan¬ 
tum correction—vacuum energy. Concerning the fundamental scalar held, it faces the same 
problem of “naturalness”. In Minkowski space time, vacuum energy can be safely ignored 
by choosing normal ordering in quantum held theory. In this case, it neither inhuences the 
dynamics nor there is a known way to measure it in local experiments. However, in curved 
space time, vacuum energy adds to the energy momentum of matter and it contributes to 
the dynamics of spacetime geometry. In the important case of Freidmann-Robertson-Walker 
(FRW) cosmology, geometry is conformally equivalent to Minkowski spacetime and thus one 
might naively think that it has solved the cosmological constant problem. In this case, how¬ 
ever, one is left with a scalar held non-minimally coupled to matter in hat spacetime [16, 17]. 
Hence, the cosmological problem translates into an equivalent problem of naturalness of the 
scalar held. 

Leaving the aforementioned problem aside, if we opt for quintessence then a specihc 
behavior of the scalar-held dynamics is required, in order to keep the thermal history of the 
universe intact. Moreover, if quintessence has originated from inhation in the early universe 
in an attempt to attribute both accelerating phases in the same scalar held (quintessential 
inhation) [18-22], then the scaling behavior of the scalar held is required in order to comply 
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with the nucleosynthesis constraint. The exit from the scaling regime to the late-time 
acceleration can be caused by invoking the non-minimal coupling with massive neutrino 
matter [23-25]. It is really mysterious that the neutrino mass is around the mass scale 
associated with dark energy. Once the non-minimal coupling to neutrino matter is invoked, 
the mass of the neutrino gets connected to the minimum of the effective potential and 
thereby to dark energy [22]. The neutrino mass modihes the cosmological evolution at both 
background and perturbation levels [26, 27]. Amongst others, the neutrino mass shifts the 
time of matter-radiation equality, and the free-streaming massive neutrino suppresses the 
growth of matter density perturbation within the free-streaming length scale [28, 29]. It is 
well known that particle-physics experiments give rise to constraints on neutrino masses. 
On the other hand, cosmology provides an independent probe, which is however sensitive to 
the new degree(s) of freedom over and above the standard model of particle physics. The 
massive neutrino matter crucially affects the matter power spectrum, and thus it can give 
rise to constraints on neutrino masses. In particular, the matter-density power spectrum is 
damped on small scales by massive neutrinos. 

Having in mind the above discussion, it is both interesting and necessary to investigate 
the observational constraints on varying neutrino-mass cosmology. Although observational 
constraints analysis on varying-mass neutrino cosmology has been performed in the past (see 
for instance [30-35]) the crucial new feature of our investigation is the incorporation of the 
non-minimal coupling. In particular, we shall investigate quintessence models based upon 
the tracker held with massive neutrino matter minimally as well as non-minimally coupled to 
the held. Moreover, another novel feature is the use of the recently released Planck and SDSS 
data. The plan of the manuscript is the following: In section II, we consider dark energy 
scenarios based upon scalar held models, with inverse power law and double exponential 
potentials with minimally coupled neutrino matter. Section III is devoted to the dynamics 
of a scaling held coupled to massive neutrino matter. In section IV we present our results 
on the observational constraints on the above models. Finally, we give our conclusions in 
section V. 
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II. QUINTESSENCE MINIMALLY COUPLED TO MATTER 


The quintessence scenario is one of the main paradigms for the description of the late¬ 
time acceleration [36]. In this set up, along with cold dark matter and baryonic component, 
we shall also be interested in considering massive neutrino matter. As mentioned above, 
this class of models can give rise to quintessential inflation, such that the dark energy scale 
emerges naturally in the scenario. 

Before proceeding to non-minimal coupling with neutrino matter, we first briefly consider 
standard quintessence with minimally coupled massive neutrino matter. In order to set the 
notations, let us consider the following standard action 


iS = y 


/Lf2 1 


+ Sm + Sr , 


( 1 ) 


where Sm, r denotes the action of ordinary matter including non-relativistic and relativistic 
one. Additionally, in order to describe the background cosmological evolution, we consider 
the spatially flat FRW geometry, in which case the action (1) leads to the following evolution 
equations 


3i7^Mp[ — Pm + Pr + 5 

(2) 

( 2 H ilLpj = — (pm + Pr + P 0 ) 5 

(3) 

• dV 

(4) 

0 -|- 3770 H—j— — 0 , 
d(p 


where p^j, and denote the energy density and pressure of the quintessence field. In this 
case as usual, the dark energy sector is attributed to 0 and the dark-energy equation-of-state 
parameter reads w^e = = p^/p^. 

As mentioned in the introduction, specific features of scalar field dynamics such as tracker 
solutions are of interest in cosmology. In this case, once the conditions at the present epoch 
are set properly, cosmic evolution is not sensitive to the initial conditions. In what follows, 
we shall try to alleviate the problem associated with the said choice, for model parameters 
attributed to tracker solutions. 

To make the point, we shall consider tracker solutions in the models with inverse power- 
law and double-exponential type potentials, though the underlying features are common to 
any tracker model. Let us begin with the inverse power-law type potential [13, 37-39] 
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where Vq is a mass dimension-four constant. In this case, the slope of the potential is given 
by = —MpiV'/V = n/(j), which is large at early times and diminishes at late times when 0 
acquires large values. Consequently, the held might mimic the background {scaling solution) 
at early epochs and could exit to slow roll at late times, giving rise to de Sitter attractor a 
la a tracker solution. As explained in the Appendix, it is difhcult to achieve tracker in this 
case for generic values of model parameters. For instance, for Vq ~ Mp, one requires much 
larger value of n (see Ref. [20] for details). It turns out that it is much easier to realize 
tracker in a model with the double exponential potential. 

Let us consider the potential of the form [40-42] 

R(0) = Vo , (6) 

where Vq, Ci and C2 are constants. For (/)j 3> 0 at the initial time (A^ = Ina —)■ — oo), the 
second term in (6) can give rise to scaling regime for generic values of C 2 > 9 consistent with 
nucleosynthesis constraints of = 3(1 -|- Wb)/cl < 0.045 [33, 43], where Wb is the equation- 
of-state of the background huid, i.e., Wb = 1/3 and 0 in radiation and matter dominated 
epochs, respectively. On the other hand, for late times, when 0 approaches the origin, (6) 
exhibits a minimum with Vmin ~ Vo for C 2 ~ ci. The latter can suit our requirement if we 
choose Vq ~ Per- 

In order to present the aforementioned behavior in a more transparent way, we numer¬ 
ically evolve the evolution equations and we depict the corresponding evolution in Fig. 1. 
Indeed, in the upper graph of Fig. 1, we depict the corresponding evolutions for the energy 
densities of radiation {pr), matter (p^) and quintessence held (p^), normalized to the matter 
energy density pm at present, as functions of A^ = Ina, for three choices of Vq, ci and C 2 . 
In the lower graph of Fig. 1, we depict the corresponding evolution for the quintessence 
equation-of-state parameter Because the cosmic evolution is insensitive to the initial 
conditions in the exponential potential, we can choose initial values in a broad parameter 
space. In the calculation, 0/Mpi = 30/c2 and 0 = 0 are set at A^ = —20. As we observe, 
the density ratio p^/Pm , depending on C 2 in the early universe, does not change in radiation 
and matter dominated eras. Finally, the behavior of 0-held at early times, can be estimated 
by combining (4) and (6), leading to = Wb- Hence, as the scale of p^p is close to Vq, the 
quintessence held evolves into the ACDM-like stage, where wp, —1. Note that because 
the late-time behavior of the double exponential potential is almost the same as that of the 
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FIG. 1: Left: Evolutions of various energy densities, normalized to matter energy density 
at present (oq = Ij, as functions of N = In a, for the minimally coupled scenario with the 
double-exponential potential V = pr is the gray-dashed curve, pm is the 

gray-solid curve, and p^ is plotted for three parameter choices, namely (ci,C 2 , Vo/p® )= (4,50,2) 
(blue-solid), /ci, C 2 , Vo/p®/= (4, 50, 20) (green-dashed) and (ci,C 2 ,Vq/ pm^ )= {4,10,2) (red-dotted). 
Right: The corresponding evolution of the quintessence equation-of-state parameter w^. We have 
used p^'^ /p^ = 3 X 10“^ as the boundary condition. 

TABLE I: List of priors for parameters and allowed regions with 95% C.L., and E((/>) = 
Vb ^ 20. 


Parameter 

Prior 

Our result (95% C.L.) 

ACDM (95% C.L.) 

Baryon density 0.5 < 100Lli,h? < 10 

lOOLlbh"^ = 2.201[];[]5 

lOOClbh 

2 9 99+0.04 

. — Z.ZZ_Q Qg 

CDM density 

10“^ < LLch^ < 0.99 

Q.ch^ = 0.118 ±0.003 

Llch^ -- 

-Oil 7+0-004 

- U.ii r_Q QQ2 

Neutrino mass 

0 < < 1 eV 

Smp < 0.200 eV 

Umiy 

< 0.198 eV 

Spectral index 

0.9 < Us < 1.1 

Us = 0.961 ±0.011 

Us = 

0.963l°+? 

Potential 

10 < C2 < 40 

C2 > 24.6 


- 


ACDM model, this scenario is hardly distinguished from ACDM. We perform the CosmoMC 
package [44, 45] to extract the observational constraints and present the results of these two 
scenarios with one massive neutrino along with the other two being massless in Fig. 2 and 
Table. I. The obervational dataset will be introduced in Sec. IV. 

It is clear from the above that the power-law potentials inevitably lead to a dark-energy 
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FIG. 2: Likelihood contours of the mass sum of the three neutrino species Timy and energy density 
ratio klyh? versus the CDM physical density Qch'^, for the minimally coupled scenario with the 
effective potential V{(f>) = Vq where the inner and outer curves represent 1 

and 2a confidence levels, respectively. 

dominated epoch, and this (/)-dominated era occurs when is of the same order as pm- 
However, for generic values of model parameters, the scenario is severely constrained by 
observations. On the other hand, the double exponential potential leads to a different 
behavior, in which a natural evolution in the radiation, matter and late-time acceleration 
regimes can be realized, provided that we choose appropriately the model parameters. In 
particular, the appearance of the dark energy dominated era crucially depends upon the 
scale of the potential Vq, which should be chosen around pcr- The latter brings back the 
original problem associated with the cosmological constant. 

In this section, we have performed the analysis in tracker models that have generic features 
related to the insensitivity of evolution to initial conditions. In this case, our results are 
close to the ACDM paradigm [9]. This is not surprising since we are dealing with the slowly- 
rolling scalar held, which mimics the cosmological constant at late times. Therefore, we do 
not have signihcant improvements compared to the predictions of the ACDM scenario. It 
is clear that a successful model of quintessence should maintain the scaling behavior, with 
an additional mechanism of exit to dark energy that could avoid the above disadvantage. 
This mechanism is provided by the non-minimal coupling between quintessence and massive 
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neutrino matter. In the section to follow, we construct and analyze a model with the above 
features. 


III. NON-MINIMALLY COUPLED MASSIVE NEUTRINO MATTER AND 
GRACEFUL EXIT TO DARK ENERGY 


As mentioned above, one should look for a scenario in which the scale of dark energy 
emerges somewhat naturally. In this section we briefly review a model in which a key 
mechanism is the nonminimal coupling of the scalar field with varying-mass neutrinos [22]. 
In this case, Vq gets linked to the energy density of massive neutrino matter. 

The relevant action in the Einstein frame is written as 


S = f d‘^Xy/—g 


- + \d^ct)d^ct) + U(0) 


+*5111 + iSr -|- {C'^gap, ^1/) • 


(7) 


In the above expression, we have introduced a scalar field 0 with V (0) its potential, as 
well as matter and radiation sectors corresponding to perfect fluids. Note that due to the 
nonminimal curvature-neutrino coupling in the Jordan frame, we obtain a varying-mass 
neutrino sector in the Einstein frame, quantified by the coefficient [22] 


q2 _ ^2g2/30/Mpi 


( 8 ) 


Hence, varying the action (7) with respect to the metric and specializing to the FRW 
case, we obtain the evolution equations 


— - 0 ^ -|- V ( 0 ) -|- Pm + Pr + Pi/ , 

2H + MJ, = -ip + v(il>) -Pm-P,- 


(9) 

( 10 ) 


where pi are the energy densities of the corresponding sectors, and p* their pressures. Ad¬ 
ditionally, varying the action (7) with respect to the scalar field, we derive its equation of 
motion, which reads 

■ dV B 

+ 3770 -|- —— = — 3pi/), (11) 

d0 A7pi 

and thus the evolution equation for the neutrinos becomes 


py + 3i7(p^ + p^) = - 3p^). 

Mpi 


(12) 






Moreover, concerning radiation and matter sectors, we have the standard evolution. Addi¬ 
tionally, for the neutrino sector we have 

, (13) 

and thus the neutrino pressure behaves as radiation during the early times, while it 
behaves as non-relativistic matter during late times [24] where non-minimal coupling builds 
up to induce non-zero neutrino masses. In this case, at late times, one obtains an effective 
dark energy sector, consisting of the scalar field as well as of the varying-mass neutrinos. In 
this case, the energy density and pressure of dark energy are given by 


Pde = Pu + P<t> — Pu + 2^^ + ^(0)) 

(14) 

o 

H 

III 

+ 

+ 

toi 

to 

1 

(15) 


which obey the standard continuity equation. Finally, as usual it proves to be convenient to 
introduce the dimensionless density parameters for radiation, matter, neutrinos and scalar 
field as 

> (i = m,r,z/,0) (16) 

respectively, and therefore from (14) we obtain 

UdE = -|- fliy . (17) 

In order to obtain a viable model, we have to consider a specihc ansatz for the potential 
V (0) which is dictated by the requirement of nucleosynthesis. Hence, in order to illustrate 
the neutrino matter effect we consider the scaling potential, namely 

1/(0) = . (18) 

The exponential potential in the non-minimally coupled quintessence model has been inves¬ 
tigated with a < \/2 [32, 35], satisfying the slow-roll potential condition. On the other hand, 
the scaling potential, a S> 2, has also been discussed to address the late-time acceleration 
occurring due to the effective potential of the non-minimally coupled neutrino matter [33]. 
Our work follows the later “naturalness” scaling potential with current observations. Thus, 
the scalar held equation of motion (11) becomes 

0 -|- 3i70 -|- Veff^(j) = 0 , (19) 
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where 


is the derivative of the effective potential in which the scalar held moves, constituted of the 
original potential as well as of extra terms contributed by massive neutrino matter coupled to 
scalar held. The above cosmological scenario proves to exhibit a very interesting behavior, 
since this ehective potential develops a minimum. We note that the latter provides an 
interesting mechanism for the unihed description of inhation and late-time acceleration [24, 
46]. Indeed, we easily check that the minimization of 14// requires, 

^ Mpi ^ f Vo ^ 

(Pmin = -In - ; P = a'f (21) 

which leads to 

y™" = (l + 7)p.(0m.n), (22) 

where we have ignored the neutrino pressure at late times. We observe two features from 
(22). Firstly, it is clear that 147/” very sensitive to 14 and that the minimum of 

the ehective potential is mainly dehned by the numerical value of neutrino matter density 
at the present epoch. Secondly, the late-time dark energy density depends on the model 
parameters 7 . It is interesting to note that u^e = —1 —2if/3Fr^ = — 7 /(l-|- 7 ) at the present 
epoch. There is no much hne tuning in the numerical value of 7 and dark energy density in 
this case gets naturally connected to massive neutrino matter density at the present epoch. 
Note that we had introduced 7 for simplifying the analytical estimates though we do not 
need to use it any further. In our numerical elaboration we only consider a to be the free 
parameter, since it controls in matter and radiation dominated eras. 

In order to examine the above behavior in more detail, we perform a numerical elaboration 
and we depict the corresponding evolution of the various energy densities, the quintessence 
equation-of-state parameter and the neutrino masses in Fig. 3. We should emphasize that 
thanks to the tracker behavior in the high redshift regime and the emergence of the dark- 
energy dominated era at late times, the scenario can successfully be embedded in the frame¬ 
work of quintessential inflation (see Ref. [22, 24] for details). Additionally, since the early- 
time and late-time evolutions are insensitive to the initial conditions and 14 , respectively, we 
choose (4/Afpi, 4/Mpi) = (—60/a, 0) at = —20 and fix I 4 for convenience as I 4 = TpmVs, 
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FIG. 3: Upper left: Evolutions of the various energy densities, normalized with the matter energy 
density at present (oq = Ij, as functions of N = Ina, for the non-minimally eoupled scenario, 
with the effective potential /Mpi + /3{py — 3pi,)/Mpi, where pr is the gray- 

dashed curve. Pm is the gray-solid curve, and p^ is plotted for three parameter choices, namely 
a= 30 (blue-solid), 25 (green-dashed) and 20 (red-dotted), with Vq/ p^m = 7/3- Upper right: The 
corresponding evolution of the quintessence equation-of-state parameter w^j,. Lower: Evolutions of 
the neutrino masses relative to the current masses. We have used = 0.2 eV, Llmh'^ = 0.118 
and pf'"’ / Pm = 2 .58 X 10 as the boundary conditions. 


which corresponds to the dark energy scale. For this Vq choice, the hnal scalar approaches 
to the origin of coordinate, i.e., 0 ~ 0 at .2 = 0. Furthermore, we should also mention that 
in this work we focus on the simplest varying neutrino-mass model, with only one neutrino 
flavor having the held varying mass, while the other two are treated as massless. 

As the matter density perturbation is suppressed by the free-streaming massive neutrino, 
the neutrino mass inhuences the matter power spectrum. However, since small neutrino 
masses are preferred for non-minimally coupled neutrino matter in the early time (see also 
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the lower panel in Fig. 3), the suppression in the matter power spectrum due to massive 
neutrinos is minimized. In the upper panel of Fig. 4, we depict the matter power spectrum 
with EttIj, = 0.04 and 0.15 eV, where the measured data points are from the Sloan Digital 
Sky Survey (SDSS DR7). In the lower panel of Fig. 4, we plot the matter power spectrum 
deviation between = 0.15 eV and 0.04 eV in coupled neutrino matter (black solid 
line) and ACDM (gray solid line) models. Compared to ACDM, we can see that in the 
nonminimally coupled neutrino matter model the suppression of the matter power spectrum 
is minimized due to the free-streaming neutrino. 

As shown in Refs. [26, 27], the scalar 0 within a large scale neutrino lump largely deviates 
from that at the background level 0 , leading the neutrino mass to become negligible inside 
the neutrino lump, and thus the perturbation to become non-linear. In this work we desire to 
avoid such a non-linear region and we assume that the neutrino mass rrii, and the scalar 0 are 
both homogeneous, and therefore the perturbation of neutrinos in the nonminimally-coupled 
scenario behaves as that in the ACDM paradigm: 

s)‘ 

Finally, since in the scenario at hand the dark energy sector is attributed to the 
combination of the quintessence held and the varying-mass neutrinos according to (14) 
and (15), for completeness in Fig. 5 we depict the evolution of along with 

wde = (P 0 + Pu)/{P 4 ) + Pu) and = Pyjpy. As expected according to our previous discus¬ 
sion, at late times determines wde almost completely. In addition, from (20) we estimate 
that 

Q crR 32 q: 

H rs> _ _ 

py Eniy/eV ’ 

where we have used p^ ~ V{(j)rnin), = 0.34 and VLyh? = Emj,/94.1 eV. Combined 

with (13), the small huctuation of 0 results in a signihcant change in neutrino masses. As 
shown in Fig. 5, neutrino masses along with the equation-of-state parameters oscillate in 
—2 < Ina < —0.5, i.e. 0.7 z < 6.4, and become purely massive after 2 ; < 0.7. 

As the scenario at hand leads to a successful description of the universe history at the 
background level, it is necessary to confront it with observations-related matter perturba¬ 
tions, as the varying-mass neutrino could in principle lead to strong constraints in this case. 
In the following section we perform such an observational analysis in detail. 



5y = 3H 


Wy 


(1 +Wy) [9 


(23) 
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FIG. 4; Left: Matter power spectrum of the non-minimally coupled neutrino matter model with 
^eff,cj> = —oiVoe~°‘^^^P^/Mpi-\-f3{py—3py)/Mpi, where Yim^ = 0.04 eV (solid) and 0.15 eF (dashed), 
Vo/P^m = 7/3 os the boundary condition and (a,Tjmy/&Y,Q.ch?)= (25,0.2,0.118,2.20). 
Right: the deviation of the matter power spectrum AP^ = [Pqi^ — 7 ^ 004 ) 1^0 04 , coupled neutrino 
matter model (black solid line) and AP^^^^ = in ACDM model 

(gray solid line), where Po.i 5 ,o .04 = P{P) with the neutrino mass = 0.15 and 0.04 eF. 

IV. OBSERVATIONAL CONSTRAINTS 


We use the program CosmoMC [44, 45] to extract the observational constraints on the 
scenario of the previous section. An important point is that one should go beyond the 
background evolution, since the massive neutrino effects will also arise from the mat¬ 
ter density perturbation. Therefore, in our analysis we incorporate the matter power- 
spectrum data sets, which include cosmic microwave background (CMB) from Planck [9] and 
WMAP [47], baryon acoustic oscillation (BAO) from Baryon Oscillation Spectroscopic Sur¬ 
vey (BOSS) [48, 49], Type-la supernova (SNIa) from Supernova Legacy Survey (SNLS) [50], 
and matter power spectrum from Sloan Digital Sky Survey (SDSS DR4) [51] and Wig- 
gleZ Dark Energy Survey [52, 53]. The details of the htting procedure can be found in 
Refs. [44, 45]. In Fig. 6, we depict the 2D likelihood contours for the mass sum of the three 
active neutrinos at present {z = 0) and the model parameter a, versus the physical density 
of cold dark matter (CDM) The scenario at hand is in agreement with observations, 

and indeed quantities like the present dark-matter and dark-energy density parameters have 
similar ranges as those in ACDM cosmology [9]. However, the neutrino mass sum in the 
model is enhanced to 0.28 eV and the allowed window is also signihcantly relaxed such that 
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FIG. 5: Evolutions of the equation-of-state parameters of the total dark-energy sector wde 
( blue-solid), of the quintessence field (green-dashed), and of the neutrinos Wy (red-dotted), 
as functions of N = In a, for the non-minimally coupled scenario, with the effective potential 
^eff, 4 > = —cxVoe~°^^^^P^/Mpi + I3{pu — 2>pu)/Mpi, where we have used a= 25 and Vq/ p^m = 7/3 o-s 
well as TiItIu = 0.2 eV, Qmh'^ = 0.118 and p^^/Pm^ = 2.58 x 10“^ as the boundary conditions . 



FIG. 6: Likelihood contours of the sum of the mass of the three neutrino species in eV 

(upper) and the potential parameter a (lower) versus the CDM physical density Llff?, for the 
non-minimally coupled scenario with the effective potential 14 //,,/, = /Mpi + j3{piy — 

'ipu)/Mpi, where the inner and outer curves represent 1 and 2a confidence levels, respectively. 


14 






















TABLE II: List of priors for parameters and allowed region with 95% C.L., and Veff^tj, = 
+ /3ip, - 3p,)/Mpi. 


Parameter 

Prior 

Our result (95% C.L.) 

ACDM (95% C.L.) 

Baryon density 

0.5 < < 10 

lOODfe/i^ = 2.19l|];[|^ 

lOODft/i^ = 2.221°;°^ 

CDM density 

10-3 < < 0.99 

= 0.119 ±0.003 

Q ^2 _ Q 1 i'v+0.004 

— U.il/_QQQ2 

Neutrino mass 

0.03 < Tirrii, < 1 eV 

= 0.278toiil eV 

Srui, < 0.198 eV 

Spectral index 

0.9 < Us < 1.1 

Us = 0.963l°;°J? 

n, = 0.963;°;°J? 

Potential 

20 < a < 30 

a > 22.1 

- 


SttIj, < 0.52 eV. The potential parameter a controls the ratio p<^/pc in the high redshift 
regime, which at the 2 — a confidence level is constrained to be a > 22.1, which naturally 
yields the tracker behavior so that the coincident problem could be understood. Our results 
for the non-minimally coupled scenario are summarized in Table. 11. Finally, it is worth to 
mention that the coupling 13 is tuned to control in our program, and its value is roughly 
inversely proportional to a and Smj,. In our numerical computation, we take (3 ~ SOO/Smj, 
with a = 25. This can be estimated by the relation of Eq. (24). 

V. CONCLUSIONS 

We have examined the tracker quintessence models with inverse power-law and double¬ 
exponential potentials. In the first case, it is difficult to obtain a viable scenario for generic 
values of model parameters, whereas in the second case we are forced to make a natural choice 
for the scale of the potential which brings back the problem that cosmological constant is 
plagued with. With a hope to alleviate the problem, we have considered massive neutrino 
matter non-minimally coupled to gravity. During radiation and early matter era, massive 
neutrinos exhibit relativistic behavior which implies vanishing coupling to scalar held. Only 
at late times, when massive neutrinos turn non-relativistic, their direct coupling to held 
builds up, leading to appearance of a minimum in the ehective potential of the held. In this 
picture, the minimum of the ehective potential is insensitive to the scale of the potential Vq 
and is rather given by the numerical value of around the present epoch. At late times, 
when the scalar held rolls around the minimum of the ehective potential, the dynamical 
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system enters into scaling regime, which is an attractor. Thanks to non-minimal coupling, 
the scaling solution is accelerating. Clearly, the model under consideration has an edge over 
standard quintessence with flat potential. 

We have found that the neutrino mass grows in time and its mass in the past is smaller 
than the current stage. In the early universe, the free-streaming length of the neutrino 
around the horizon scale is the same as that of the massless one. As shown in Fig. 5, the 
neutrinos become massive only at the very recent epoch. As a result, our scenario can relax 
the suppression of the matter power spectrum from the free-streaming massive neutrino 
in the early universe. The constraint on the neutrino masses is reduced in comparison to 
ACDM. In particular, in the scenario under consideration, we found that the best-fit value 
of the neutrino mass sum are around 0.28 eV, along with the allowed window extended to 
0.52 eV, which is consistent with particle physics experiments. 
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APPENDIX: THE INVERSE POWER-LAW POTENTIAL 


For completeness, let us investigate the case of the power-law potential (5). In this case, 
the evolution equations can be explicitly integrated in the background (radiation/matter) 
dominated regime. The asymptotic behavior of the solution in the radiation and matter 
dominated eras is given by a power-law behavior with 

—2 -|- nwb 



Eq. (25) implies that < Wb for n > 0 in the radiation and matter dominated epochs, i.e. 
P<p/Pm increases, allowing finally to catch up with dark-energy domination. For moderate 
values of n ~ 1, the held might take over the background and as a result pp, might obtain 
the observed value of dark energy at the present epoch, whereas has not yet reached the 
desired value. In the upper graph of Fig. 7, we depict the corresponding evolutions for the 
energy densities of radiation (p^), matter (p^) and quintessence held (p<^), normalized with 
the matter energy density p^m at present (uq = 1), as functions oi N = In a, for the case 
n = 4 with three choices of Vq. In the lower graph of Fig. 7, we show the corresponding 
evolution for the quintessence equation-of-state parameter w^p. From this plot we can see 
that there exist two regimes with constant values of wp, ~ —1/9 and Wrp ~ —1/3, in the 
radiation and matter dominated era, where the held tracks the background. Clearly, the 
equation-of-state parameter w^p does not reach the required value at the present epoch. 

The above situation can be remedied by taking large values of n, for which wp, Wb. In 
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FIG. 7: Left: Evolutions of the various energy densities, normalized with the matter energy density 
p® at present (oq = Ij, as a function of N = In a, for the minimally coupled scenario with the 
inverse power-law potential Vq ■ Pr is the gray-dashed curve, pm is the gray-solid curve, 

and Pfj, is plotted for three parameter choices, namely (n, Vq/ p^m (4,1) (blue-solid), (n,Vo/pm'^)= 
(4,10^) (green-dashed) and (7i,Vo/pm^j= (4,10^) (red-dotted). Right: The corresponding evolution 
of the quintessence equation-of-state parameter w^. We have used p^'^/p^ = 3 x 10“^ as the 
boundary condition. 

this case, the inverse power-law potentials reduce to the exponential form, leading to scaling 
solution. The late-time exit from the scaling regime is guaranteed by the shallow nature of 
inverse power-law functions, though we need to tune the model parameters appropriately. 
Firstly, irrespectively of the quintessence model, using the slow-roll parameter rj we hnd 
that 10 eV. The slow-roll parameter e tells us that this would happen at 

the present epoch, provided that c^q/M-pi » \/2n. Hence, this implies that Vq ~ n'^pcr- In 
the case where Vq ~ Mp, one requires a large value of n. On the other hand, for moderate 
values of n we need to choose Vq ~ Per (Sahni, Sami & Souradeep [20]). 
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